All electrical measurement of the density of states in (Ga,Mn)As 
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We report on electrical measurements of the effective density of states in the ferromagnetic semi- 
conductor material (Ga,Mn)As. By analyzing the conductivity correction due to enhanced electron- 
electron interaction the electrical diffusion constant was extracted for (Ga,Mn)As samples of differ- 
ent dimensionality. Using the Einstein relation allows to deduce the effective density of states of 
(Ga,Mn)As at the Fermi energy. 



PACS numbers: 71.20.-b, 75.50.Pp, 73.23.-b 

The ferromagnetic semiconductor (Ga,Mn)As [T has 
been studied intensely over the last decade and has be- 
come a model system for future spintronics applications 
[21 [3] . With typical Mn-concentrations between 1 % and 
15 % maximum Curie temperatures of up to ^ 180 K 
have been reported [H |5] . Mn atoms on Ga-sites provide 
both holes and magnetic moments. For Mn concentra- 
tions larger than 1 % the impurity wavefunctions at the 
Fermi energy overlap and a metallic state forms. The fer- 
romagnetic order between the magnetic moments of the 
Mn-ions is mediated by the delocalized holes [6^ . A topic 
of current debate is whether the holes reside in an im- 
purity band, detached and above the valence band or in 
the valence band [7j. A mean field picture based on the 
latter scenario allowed to predict, e.g. Curie temperature 
[6] or magnetocrystalline anisotropics [8] in (Ga,Mn)As 
correctly. On the other hand optical absorption exper- 
iments, carried out, e.g. in Ref. [H [TO], suggest that 
even for high manganese concentrations of up to 7 % the 
Fermi energy stays in an impurity band, detached from 
the valence band, with a high effective hole mass of order 
ten free electron masses rUe [9]. However, there is also 
indication that the impurity band and the valence band 
have completely merged as discussed in Ref. ^ and ref- 
erences therein. In the present letter we make use of the 
well known quantum mechanical conductivity correction 
due to electron-electron interaction (EEI) to extract the 
diffusion constant and hence the density of states at the 
Fermi energy, N(Ep). The electrically measured values 
of N{Ep) will be compared with recent theoretical cal- 
culations. 

In ferromagnetic (Ga,Mn)As the conductivity is de- 
creasing with decreasing temperature below 10 K. This 
conductivity decrease can be explained by enhanced 
electron-electron interaction [11] . The effect of EEI arises 
from a modified screening of the Coulomb-potential due 
to the carriers' diffusive motion and depends on the di- 
mensionality of the conductor [K]. As the conductiv- 
ity decrease due to enhanced electron-electron interac- 
tion is depending on the electrical diffusion constant D, 
a detailed analysis of the conductivity decrease, different 
for different dimensionality, provides experimental access 
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TABLE I: Length I, width w, thickness t and number of lines 
parallel N of the samples. Curie temperature Tc and carrier 
concentration p were taken on reference samples from the 
corresponding wafers. Annealed samples are marked by "A". 



to the diffusion constant. Using the Einstein relation 
a — N{Ep)De'^ , with the conductivity a, the effective 
density of states at Fermi's energy, N{Ef), can be deter- 
mined. 

To investigate electron-electron interaction in quasi 
ID, 2D and 3D systems we fabricated Hall-bar mesas (2D 
and 3D) and wire arrays (ID and crossover regime from 
ID to 2D) out of several wafers, having a (Gai_a;,Mn2.)As 
layer on top of semi-insulating GaAs. The nominal Mn- 
concentration x was approx. 4 % (sample I3D and 23£)) 
and ~ 6 % (other samples). The relevant parameters of 
the samples are listed in table |lj The dimensionality for 
EEI is defined by the number of spatial dimensions larger 
than the thermal diffusion length Lt = ^JhD /ksT . In 
(Ga,Mn)As Lt « 120 - 200 nm at 20 mK, depending 
on the exact value of the diffusion constant D. Hence 
the thick Hall-bar mesas (150 nm and 300 nm) can be 
considered as quasi 3D, while the thin Hall-bar mesa (42 
nm) is quasi 2D, at least below ^ 500 mK. The small- 
est wires (42 nm and 35 nm) behave quasi ID and the 
wider wires (62nm to 242 nm) are in the crossover regime 
from ID to 2D, as is shown below. Arrays of wires with 
N wires in parallel were fabricated to suppress univer- 
sal conductance fluctuations by ensemble averaging. The 
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FIG. 1: a) Electron micrograph of a line array having 25 lines 
in parallel (Sample 2id~2dA). The width of the lines is 92 nm, 
the length is 10 /im. b), c) and d): Conductivity of the quasi 
ID line-array Im (b), the quasi 2D Hall-bar I2D (c) and the 
quasi 3D Hall-bar 230 (d) plotted versus temperature. The 
straight lines are guide for the eyes. The slope of the lines are 
given. 



Hall-bars were fabricated using optical lithography and 
wet chemical etching. For fabricating the wire arrays we 
used electron-beam lithography and chemical dry etch- 
ing. The contact pads to the devices were made by ther- 
mal evaporation of Au and lift-off. The measurements of 
the conductivity were performed in a top-loading dilution 
refrigerator using standard four-probe lock-in technique. 
To avoid heating of the charge carriers small measur- 
ing currents (25 pA to 4 nA, depending on the sample's 
resistance) and proper shielding were crucial. For each 
sample the measuring current was kept fixed for all tem- 
peratures. To suppress conductivity contributions due to 
weak localization we applied a perpendicular magnetic 
field of i? = 3 T. At B = 3 T no weak localization can 
be observed in (Ga,Mn)As dSHH] even at 20 mK. 

According to Lee and Ramakrishnan [T^] the temper- 
ature dependency of the conductivity correction due to 
EEI is depending on the dimensionality of the sample 
with respect to Lt- For ID-systems the expected tem- 
perature dependency is oc —l/VT, for 2D oc logiQ(T/To) 
and for 3D c>c ^/T. Corresponding data for ID, 2D and 
3D (Ga,Mn)As-samples, shown in figure lb)-d), confirm 
the expected temperature dependency below 1 K. Hence 
the decreasing conductance with decreasing temperature 
in (Ga,Mn)As can be attributed to EEI. 

The size of the conductivity correction due to electron- 
electron interaction is depending on the diffusion con- 
stant D in ID-systems [T^ : 
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As the conductivity correction due to EEI is also 
depending on the screening parameters p^d,2D,3D ^j^g 
needs to know the corresponding parameter to extract 
D from the conductivity correction. Only in quasi 2D- 
systems the conductivity correction is independent on the 
diffusion constant [T^: 
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Hence in the 2D case F^^ can be directly extracted from 
experiment. As already shown in previous work |llj the 
screening parameter F^^ in (Ga,Mn)As ranges from 1.8 
to 2.6 and is in excellent agreement with the screening 
parameter in Co, Co/Pt-multilayers and Permalloy. In 
these ferromagnetic metals F^^ is between 2.0 and 2.6 
[TSl [TBI El CH] • Thus using the well known parameters 
F^^ of other ferromagnetic metals is a good approxima- 
tion for F^^ in (Ga,Mn)As. In Ni and Py nanowires 
F^^ is 0.83 and 0.77 respectively [HI US]- Consequently 
F^^ — 0.80 ± 12 should be a good approximation for 
the screening parameter of quasi ID (Ga,Mn)As sam- 
ples. With this F^^ parameter we can calculate the 
diffusion constant of sample Iid using equation (1): 
D — 10 ± 3 • 10~^ m^/s. Using the Einstein relation 
this value corresponds to an effective density of states 
N{Ef) = 1.1 ± 0.3 • 10-*^ /Jm^ at the Fermi energy. In 
figure 3 N{Ep) is plotted versus the carrier concentration 
(green squares) for all investigated ID-samples. Here, 
the uncertainty in determining the carrier concentration 
is ~10 %. 

It is more difficult to estimate the value of F'^^ as no 
data are available for 3D ferromagnetic metals. There- 
fore we have to rely on theoretical predictions for the 
screening parameter: F^^ = 1.2 [2D]. Though the calcu- 
lations of the 2D screening parameter (F^^ = 2.3 [5U]) 
agree well with the experimental values of different fer- 
romagnets {F"^^ — 1.8. ..2. 6), they are less accurate for 
ID systems. For ID systems F^^ was calculated to be 
1.6, while the typical experimental values of F^^ « 0.8 
are by a factor of 2 smaller [T5J [12] . 

Hence by using the theoretical value for F^^ we need to 
take into account an uncertainty of order 100 %. Using 
equation (2) and the theoretical value for F^^ = 1.2 
fIU\ we arrive at D = 2.2 • 10~^ m^/s for sample 23X). 
Using the Einstein relation this corresponds to N{Ef) = 
2.0 • lO''® /Jm"^, with a high uncertainty of approx. 200 
%, as -D is depending quadratically on F^^ . 

The different temperature dependence of EEI in ID 
(Act cx -1/VT) and 2D (Act cx InT) together with the 
temperature dependence of the thermal length Lt = 
y/hD/ksT allows to use another scheme to extract the 
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FIG. 2: a) Conductivity change from 1 K to 22 mK of four 
line arrays (sample lii3_2DA...4i_D_2D A) plotted versus the 
wire width. The red line is the best fit of the data to equa- 
tion Q . The blue and green lines are fits using equation ffl 
and a diffusion constant of 12 • 10~^ m^/s and 6 • 10~^ m "s 
respectively. In the inset the conductivity change of four wire 
arrays (sample liD-2DA...4iD-2i3A) is plotted versus l/\/r. 
The slopes are guide for the eyes, b) Conductivity change 
of the four line arrays (sample 1i_d_2I3A...4ii3_2I3 A) plotted 
versus temperature. The red lines are calculated using equa- 
tion (4) and the parameters obtained by fitting the data in 
a). 



diffusion constant and hence N(Ep). By measuring the 
dimensional crossover, i.e. the change of the temperature 
dependence of the conductivity correction as a function 
of the sample size one can fit both, diffusion constant 
and screening parameter independently. Here we used 
the crossover from ID to 2D to determine D. For this 
experiment wire arrays with wire widths, ranging from 62 
nm to 242 nm (sample Ii£)-2nA...4iD^2DA in table 1), 
were patterned on the same wafer. These four wires are 
in the crossover regime between ID and 2D. This is most 
clearly seen in the inset of figure 2a and in figure 2b. In 
Fig. 2b the logarithmical temperature dependency, ex- 
pected for 2D EEI, is only describing the widest wire at 



high temperatures satisfactorily, while the 1 / Vt depen- 
dency, expected for ID EEI, only holds for the smallest 
wire at low temperatures in the inset of Fig. 2a. In the 
crossover regime the conductivity correction due to EEI 
is given by an interpolation formula f21J: 
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with one screening parameter F and To, the lowest tem- 
perature. Figure 2a shows the conductivity change from 
1 K to 22 mK of all four wire arrays. The conductivity 
change increases markedly with decreasing wire width. 
To extract the characteristic parameters we fitted the 
data using equation (4) with D and F as free parame- 
ters. The diffusion constant affects essentially the width 
dependence (x-scale) while the screening parameter F 
shifts the curve on the y-scale. Hence the fit is unique 
and allows to extract D and F independently. The best 
fitting result was obtained by using D — 9-10~^ m^/s and 
F ~ 1.08 (red line). To illustrate the sensitivity of the 
fitting procedure on D we also plotted equation (4) using 
D = 12- 10-5 mVs (blue line) and D = 6 • 10"^ m^/s 
(green line). Here F was the free parameter to optimize 
the fit. Both traces (blue and green) describe the experi- 
mental data less satisfying than the red trace. Hence the 
measurement of the dimensional crossover from ID to 2D 
results in 13 = 9 ± 1.5 • 10"^ m^/s. In figure 2b the con- 
ductivity change with respect to 22 mK is plotted for all 
four wire arrays (sample 1id-2D^ to 4i£)_2dA) versus 
temperature. The red lines are the calculated conduc- 
tivity correction given by equation (4). The parameters 
used were D = 9 - 10^^ m^/s and F — 1.08. The conduc- 
tivity correction in the whole temperature range from 22 
mK to 1 K of all four wire arrays is perfectly described 
by using only these two parameters D and F. Also in 
the crossover regime from ID to 2D the observed screen- 
ing parameter F = 1.08 is in excellent agreement with 
the screening parameter observed in Co {F = 0.95) [TS] , 
This underlines the universal character of the conductiv- 
ity correction due to EEI. From the obtained diffusion 
constant we can estimate the effective density of states 
using the Einstein relation: N{Ef) ^ 1.6 ± 0.3 • lO''^ 
/Jm^. 

To check the consistency of both presented meth- 
ods, we can also treat the four samples (lii:i_2DA to 
4iD-2i>A) as quasi ID at low temperatures (indicated 
by the straight lines in the inset of figure 3a) and cal- 
culate the diffusion constant using equation (1) and 
F^^ = 0.8 as described above. When doing so we obtain 
D = 8.4±2.5-10~5 m^/s. This is in good agreement with 
the value estimated by fitting the crossover from ID to 
2D and hence the two methods are consistent. 
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FIG. 3: Density of states at Fermi's energy plotted versus 
carrier concentrations. Tlie squares are calculated using the 
model of Masek et al. and the model of Tang and Flatte 
respectively. The solid line gives N{Ef) calculated using a 
6 X 6 fc ■ p model with parameters for the GaAs valence band 
|25| . The filled symbols are measured using EEI in ID (green 
triangle), 3D (blue triangle) and in the crossover regime from 
ID to 2D (blue diamond). 



model of Tang and Flatte, the relative high uncertainty of 
the experimental results does not allow to give a definite 
answer which model describes the experimental results 
better. For the case that Fermi's energy lies within a 
detached impurity band no calculations for the density 
of states are available. However assuming a parabolic 
band, one finds that N(Ep) oc m*, for a given carrier 
concentration. Hence an effective mass m* » rUe, as 
expected for a detached impurity band, leads to an effec- 
tive density of states at Fermi's energy far away from the 
measured values. 

In conclusion we have demonstrated that the effective 
density of states at the Fermi energy of (Ga,Mn)As can 
be extracted from conductivity measurements, i.e. an 
analysis of the conductivity correction due to EEI. The 
measured values of N{Ep) are consistent with a picture 
that the Fermi energy is located within the GaAs valence 
band or an impurity band, merged with the valence band. 
Our experimental finding with effective masses of Irrig 
is however at odds with a detached impurity band, where 
the effective hole mass is much larger than mg. 
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Figure 3 summarizes our findings and shows the ex- 
tracted effective density of states versus the carrier con- 
centration (filled symbols). The data are compared to 
calculations based on three different models: The solid 
line stems from a 6 x 6 fc • p model with parameters for 
the GaAs valence band, calculated by T. Dietl et al. and 
based on reference |15]. The corresponding effective hole 
mass is ~ Img in the investigated range of carrier con- 
centration. In addition we performed numerical simula- 
tions which are based on a multi-band tight-binding ap- 
proach applied to disordered bulk systems using two dif- 
ferent parameter sets. The first model was derived from 
first principles calculations for (Ga,Mn)As [32] (model 
Masek). The second one describes the Mn impurities by 
a modified on-site potential and a spin-dependent poten- 
tial at the four nearest neighbor As sites which reproduce 
the experimental binding energy of 113 meV [53] (model 
Tang). A detailed description of the method and the two 
models is given in reference |24j . Neither model exhibits 
a detached impurity band for Mn concentrations larger 
than 1%. The effective masses were estimated to lie in 
the range m* = 0.4...0.6me for the considered carrier 
concentrations with only minor quantitative differences 
between the two models. Due to the increasing number 
of holes the Fermi energy moves deeper into the valence 
band with increasing disorder. The experimental data 
for the density of states lie between the predicted values 
of the models. Although the data seem to be closer to the 
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